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Abstract 

We consider a Dirac singlet fermion as thermal dark matter for explaining the X- 
ray line in the context of a supersymmetric Higgs-portal model or a generalized 
Dirac NMSSM. The Dirac singlet fermion gets a mass splitting due to their Yukawa 
couplings to two Higgs doublets and their superpartners, Higgsinos, after electroweak 
symmetry breaking. We show that a correct relic density can be obtained from 
thermal freeze-out, due to the co-annihilation with Higgsinos for the same Yukawa 
couplings. We discuss the phenomenology of the Higgsinos in this model such as 
displaced vertices at the LHC. 
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1 Introduction 


Dark matter (DM) is a main component of matter in the Universe, confirmed by various 
observations such as galaxy rotation curves, gravitational lensing. Moreover, it is sup¬ 
ported by the measurement of Cosmic Microwave Background Radiation and Large Scale 
Structure, and so on. However, we have no clue as to the DM mass and interactions other 
than gravity. Therefore, direct detection on earth, indirect detection in the sky, and direct 
production at particle colliders have been thought to be complementary for identifying the 
nature of DM. In particular, indirect detections look for the remnants of annihilations or 
decays of DM through cosmic rays coming from galaxies and galaxy clusters. 

There has recently been a lot of interest in light DM models, after new detection of X-ray 
line coming from galaxies and galaxy clusters mainly by the XMM-Newton observatory |T|. 
There are on-going debates on the possibility of explaining the X-ray line excess with 
thermal atomic transition |3], but it is worthwhile to take it to be a signal for DM and 
study the consequences of decaying or annihilating DM models [4—T|. 

Motivated by a toy model suggested by one of us [5], we consider a concrete model 
for explaining the X-ray line with the magnetic dipole moment of a weakly interacting 
massive particle (WIMP) in the context of a generalized next-to-minimal supersymmetric 
standard model (NMSSM) with an additional Dirac singlet superfield, dubbed as Dirac 
NMSSM 18,9]. Unlike the toy model where a discrete Zo symmetry for stabilizing DM 
is broken by a small amount at the cutoff scale [5], the corresponding discrete parity, 
i.e. i?-parity, in the supersymmetric (SUSY) version is assumed to be exact. Then, a 
singlet Dirac fermion or two Majorana fermions called the singlinos, introduced in the 
Dirac NMSSM, is the DM candidate, and it gets a small mass splitting for the X-ray line 
energy at 3.55 keV due to its small Yukawa couplings to the MSSM Higgses and their 
superpartners. In this case, a tiny magnetic transition dipole moment for decaying DM 
generates the X-ray line by the small Yukawa couplings of the singlinos. We regard the 
model as a SUSY Higgs-portal in the limit that gauginos, squarks and sleptons are heavy 
enough. We also include the effects of non-decoupled gauginos on the mass splitting of 
Higgsinos or singlinos. The lightness of Higgsinos and singlinos can be ensured by a chiral 
symmetry such as Peccei-Quinn symmetry while gauginos could be relatively light due to 
R-symmetry. 

The Dirac singlet fermion can keep in thermal equilibrium with the Standard Model 
(SM) particles at freeze-out, due to the co-annihilation with the Higgsino-like fermions. 
Consequently, we show that the correct relic density can be attained, being compatible 
with the X-ray line. In the limit of heavy gauginos, the mass splitting of Higgsino states 
is about keV scale as for the singlino fermions, so neutral or charged Higgsinos decay into 
a singlino +Z*/W*, leaving a displaced vertex due to small Yukawa couplings of singlinos. 
We discuss the possibility of discovering Higgsinos at the LHC in this new topology. 

The paper is organized as follows. We begin with the model description of the SUSY 


1 


Higgs-portal for the low-energy mass spectra of DM and Higgsinos. Then, we present 
the results of the magnetic transition dipole moment between two singlinos at one loop 
in onr model and show the parameter space that is consistent with both the energy and 
flux for the X-ray line. In turn, we discuss the bound from the DM relic density and its 
compatibility with the X-ray line. Finally, conclusions are drawn. 


2 Supersymmetric Higgs-portal 

The dark sector couples to the SM particles only through the Higgs and its superpartners. 
As an example, we consider an extension of the Higgs sector in the MSSM with a Dirac 
singlet chiral superfield containing two additional singlet superfields, S and S. We assume 
that the gauginos as well as the superpartners of quarks and leptons are sufficiently heavy 
so that they are not relevant for our discussion. Meanwhile, we also discuss the effects of 
non-decoupled gauginos in this section. 

The part of the superpotential containing only Higgs doublets, H u and H d , and the 
singlet chiral superfields are 

W 0 = A s SH u H d + XgSH u H d + M s SS + \i H H u H d + fi s S + ^S. (1) 

In this model, the Dirac singlet chiral superfield communicates with the SM only through 
the Higgs and Higgsino interactions. As for the Dirac singlino, the model can also be 
called the Higgino portal. In a Peccei-Quinn (PQ) symmetric realization of the above 
superpotential, the cubic couplings for the singlet chiral superfields are forbidden, while 
the bare Higgsino and singlino mass terms and the singlet tadpole terms can be generated 
after a spontaneous breaking of the PQ symmetry by non-renormalizable interactions with 
PQ-breaking fields. 

When there is a f/(l)g global symmetry or a Z 2 symmetry distinguishing S and S, 
the operator SH u H d is forbidden. This case corresponds to the Dirac NMSSM that was 
discussed in Ref. [8j[9], where even after integrating out the singlet scalar masses with 
keeping their fermion partners, the resulting Higgs potential gets an additional quartic 
potential, \XsH u H d \ 2 , and increases the Higgs mass as compared to the MSSM. When 
the singlet symmetry is broken spontaneously or explicitly, we can write a small Yukawa 
coupling for S such that \X$\ -C A,s = 0(1). Then, the feature of the Dirac NMSSM for 
the Higgs mass can be maintained. 

On the other hand, if |As| and |A^| are comparable, the PQ symmetry only does not 
distinguish between S and S. Thus, there is no obvious reason to forbid Majorana mass 
terms such as S 2 and S 2 in the superpotential. But, if we ignore those Majorana mass 
terms under the assumption that such a flavor structure in the dark sector is determined 
by a flavor symmetry for singlinos at a high energy scale, we can explain a small mass 
splitting and a small flux required for the X-ray line for |As|, |Ag| -C 1, as will be discussed 
in the next section. 
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M 2 = 2M \, M s = 300GeV, // e ff=31 OGeV 



Mi (GeV) 


Figure 1: Masses of Higgsino-like states as a function of gaugino mass Mi. We ignored the mixing 
effect with singlinos and took tan/3 = 10 and // e ff > 0. 


The neutralino mass matrix containing the ganginos in MSSM is given in the basis 
(B, W°, Hi H°,S, 5) by 
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where v\ + = v 2 ~ (246 GeV) 2 , tan /3 = v u /vd, and the effective /i parameter is given by 

HeE = 7 h + ^s(S) + A s(S)- 

In order to keep a small mass splitting between singlinos, we take the gauginos to 
be much heavier than Higgsinos and singlinos, namely, M 12 fi e g,Ms- Then, we can 
consider only the 4x4 sub-matrix for Higgsinos and singlinos and a mass splitting of the 
Dirac singlinos is attributed to a small coupling between Higgsinos and singlinos. Then, 
keeping all the other superpartners of the SM heavy enough, we can call the model the 
SUSY Higgs-portal. 

In the limit of M\ 2 /r e ff> Ms, the mass eigenvalues for Higgsino-like neutralinos are 
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while those for singlino-like neutralinos are, for Ag, X § -C 1, 
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Consequently, the mass differences between the nearest neutralinos are 
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Ui-MJ) 2 VJWi + M 2 ,' ’ 


where 


A± = ^(As ± Ag). 


( 6 ) 

(7) 


We note that as far as Ag and Ag are comparable, Am .34 is positive so that Xa is 
the Lightest Supersymmetric Particle (LSP) and is Next-LSP in our model. When 
the singlino mass splitting is about a few keV and g e s > M s ~ 100 GeV, the Yukawa 
couplings, Ag and Ag, should be of order 10 ~ 5 and the gaugino masses should be greater 
than about 1 TeV, unless there is an accidental cancellation^] In Fig. [lj we have illustrated 
the masses of Higgsino-like neutralinos as a function of the gaugino mass. For gaugino 
masses being greater than 1 TeV and DM mass being 300 GeV, the Higgsino mass splitting 
A/n 21 is less than 6 GeV. In Fig. [2] we show the parameter space for the Yukawa couplings 
and the mass parameters satisfying the mass splitting between singlino-like neutralinos, 
|Am 34 | = 3.55keV, in blue dashed line. 


Tn the case of cancellation, the Yukawa couplings, As and Ag, 
to increase the Higgs mass 
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can be of order one so they can be used 
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Figure 2: (Left) Parameter space of // e fr — Ms vs |As|- (Right) Parameter space of /i e g — Ms vs 
1 — A 5 /A 5 . In both figures, the parameter space explaining the X-ray line is shown between two 
red solid lines and the X-ray line energy at 3.55 keV is obtained for the blue dashed line. We took 
m H ± = 1 TeV, Mi = O. 5 M 2 = 3TeV and tan/3 = 10. 


The mass eigenstates are found by H% = Nnx°, N i 2 x°i, S = Y,i N isX° and 

S = Ni4X°i- For M\ 2 > Refr, Ms, and A s , As < 1, they read 
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The chargino mass matrix in the basis (ID+, LA+, ID , H d ) is 


Mp=| D 

' /^eff 

Then, for M 2 fi e g, the mass eigenvalues for charginos are 

Meff + V 2 sin(2/5) 


( 12 ) 


/J-efF 


m~± = M 2 + 
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(13) 


The mass difference between the lighter Higgsino-like neutralino and the lighter chargino 
is 

m xt ~ m x°i = (( X + sill ( 2 / 3 )) + i 1 + (! - 2sgn(/i ef f)) sin(2/3)) . (14) 

In Fig. [lj we have also shown the masses of Higgsino-like chargino as a function of the 
gaugino mass. In this example, the mass difference between the lighter Higgsino-like neu¬ 
tralino and the Higgsino-like chargino is less than 2.5 GeV for gauginos being heavier than 
1 TeV. 


Before closing the section, we remark on the scalar sector of the SUSY Higgs-portal. 
Due to the small Yukawa couplings of the singlinos, their superpartners, singlet scalars, 
have only a small mixing with the MSSM Higgs fields so the Higgs sector is MSSM- 
like. Moreover, it would be hard to produce singlet scalars at the current LHC at a 
detectable level. On the other hand, singlet scalars may induce the self-annihilation and 
co-annihilation of DM through the s-channels. The self-annihilation is suppressed due to 
small Yukawa couplings outside the resonance, while the co-annihilation is sizable enough 
to keep DM in thermal equilibrium, as will be discussed in Sec. [4j 


3 Magnetic dipole moments and the X-ray line 


The magnetic (transition) dipole moments can be obtained for either Majorana (5,10 
or Dirac [TT] singlet DMj^jln the SUSY Higgs-portal model, the heavier singlino yji is a 
Majorana fermion that has an almost degenerate mass with the lighter singlino xt- As 
gauginos, leptons, and squarks are assumed to be decoupled in our model, only charged 


2 Similar studies on magnetic dipole moments have been done in light of the X-ray line in Ref. 6 7. 
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Higgs and VH-boson loops contribute to the magnetic dipole moment. In order to compute 
the magnetic transition dipole moment for singlinos, we choose the non-linear R^ gauge 12 , 
in which 7 — W ± — G T interactions are absent for the charged unphysical Goldstone boson 
G ± . I 11 this case, we have to deal with only the Yukawa interactions for Goldstone bosons, 
thus simplifying the calculation. 


The singlino Yukawa interactions with charged Higgs ( H ± ) and charged Goldstone ( G ± ) 


are 


-C s = sin (3 X 2 Pl(^sS + \§S)H + cos (3 xtP l(^sS + X§S)H + 

— cos (3 X 2 Pl(^sS + A sS)G + sin (3 xt Pl(^sS + A s;S)G+ + h.c. 


(15) 


Then, from Eqs. (10) and (11), we get 


-C S = X2 (/ 3 lPl + IsrPr^H- + xiiU lPl + UrPr)xIH~ + h.c. 
+ (H~ — > G~ , sin (3 —> — cos f3, cos (3 —> sin (3 ) + h.c. + • • • , 


(16) 


where 


f 3L = iA_ sin/3, 
fsR = —*A_ cos/3, 
fm = A+ sin 13, 

fm = A + cos/3. (17) 


On the other hand, the interactions between the IH-boson and singlino-like neutralinos 
come from the mixing with Higgsinos, given as follows. 

- Cv = WYPlHIW- + %YPl8°W+ + h.c. (18) 

Then, from Eqs. Q and Q, we get the singlino-like interactions to the IH-boson as 

- £v = X2Y(93 lPl + 93rPr)x° 3 W~ + xj ^{QalPl + 9 miPn)f(lW^ + h.c. H-, (19) 

where 
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( 20 ) 
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Therefore, the magnetic transition dipole moment, generated from charged Higgs, Gold- 
stone, and IF-boson loops, is given by 


Zin^o 
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We note that due to the interchange between cos (5 and sin/3, the unphysical Goldstone 
contribution is of the same magnitude but the opposite sign as compared to the charged 
Higgs contribution. For /x e fr >> Ms and m#± ~ m~ ±, the H’-boson loops tend to be 
suppressed by m\ v . But, for /j e g ~ Mg, which is necessary for the co-annihilation of DM 
as will be discussed in the next section, the W-boson loops give a dominant contribution 
to the magnetic dipole moment of DM. 


We take two singlino-like neutralinos to be lighter than Higgsino-likc neutralinos and 
almost degenerate in mass. Then, the heavier singlino can decay into the lighter one Xa 
through the transition magnetic moment or the mixing with Higgsinos. The decay modes 
are X 3 —y Xa) and xl Xa 1 j v-, where neutrinos in the latter channel is due to the off-shell 
Z-boson. The energy of the monochromatic photon coming from y° —» XaI is given by 
E 1 ~ rri^o — rri^o for o 4 S> E 1 . For rn^o — m^o | <C m^o, the decay rates of the heavier 
singlino are 


„ „ e .mrrmo / o 

r(x§ -2 y lT ) = I! x. 


7 r 
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1 v 2 cos 2/3 
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Due to an extra factor (A'/b 34 ) 2 , the decay rate for —>• x° A vv is suppressed as compared 
to the one for X3 XaT- Thus, it is sufficient to consider only the decay mode y 3 — > XaT 
to determine the decay rate of DM. 

Suppose that the heavier singlino constitutes a fraction of the total DM by r = 
D^o/Ddm- Then, for the X-ray line at 3.55 keV, we need to take the necessary value 
of the lifetime of the heavier singlino to be r ~0 = 0 . 20 - 1.8 x 10 28 sec (7.1 keV/m^o)r |TJ[5j, 
which is equivalent to = 0.36 - 3.3 x 10 -52 GeV (m^o/7.1 keV)r _1 . For comparably 
small A s and \§, and a small mass splitting between singlinos, two singlinos contribute 
to the relic density equally, that is, r = 1/2. In Fig. [2j we show the parameter space for 
the mass splitting /i eff — Ms vs A 5 or 1 — A^/Ag, satisfying the X-ray line flux (in the 
region between two black solid lines) and the X-ray line energy (in the blue dashed line). 
Therefore, the singlino Yukawa couplings of order 10 -5 required for the X-ray line energy 
is consistent with the X-ray line flux, as far as both Yukawa couplings are of similar size, 
that is, Ag/A 5 ~ 0.97 — 0.99. 


4 Dark matter relic density 


Depending on the singlino Yukawa couplings to Higgsinos, A 5 and Ag, the singlino DM 
may be in thermal equilibrium with the SM particles due to self-annihilation and/or co- 
annhiliation with charged and neutral Higgsinos |7|. The annihilation channels for singlinos 
are ^ ff,ZH°(h °), W + H~,W+W~, y°X ?, 2 ff and x*xt ff (i,j = 3,4). 

In the case with small A 5 and X§, the self-annihilation cross sections would be too 
small to make DM in thermal equilibrium with the SM particles. However, DM can keep 
in thermal equilibrium until freeze-out, through the scattering off of the SM particles or due 
to a sizable co-annihilation with neutral or charged Higgsino by crossing symmetry jl31. In 
this case, we can obtain a correct relic density for DM, after Higgsinos are decoupled from 
the SM bath and decay into DM. Therefore, we need A s, Ag > ICY 5 for thermal DM jTj . 
The mass splitting between singlinos is 3.55 keV, so it can be ignored in computing the 
relic density. 

The relic abundance is given by 




DM — 


8.8 x 10 ” 11 GeV ~ 2 
\/<h ST, dx ( a ‘«v)x~ 2 ' 


(26) 


where g* is the effective number of relativistic degrees of freedom at freeze-out and x = 
m DM/T which read Xf ~ 20 at freeze-out temperature. The effective cross section is a 
weighted average of the annihilation cross sections for the co-annihilating particles and is 
given 1 7,13 1 by 


(cr eff u) = 


E 


i,j < 7 ij w i w j 

(i>o 2 


(27) 
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^dm^ 2 (Planck 3 cr) 



^DM(GeV) 


Figure 3: Parameter space of m dm vs /r e ff — Ms, satisfying the relic density within Planck 3er 



where 


Wi = ( 1 + A i) 7 e 


3/2 —xA i 


A _ Mi ~ ^DM 

= 5 

^DM 


(28) 


and (crijv) = cr^-x n with n = 0 (1) for s-wave (p-wave) annihilation. 


For small A s and \$, the effective annihilation cross section is dominated by the anni¬ 
hilation of neutral or charged Higgsino so we find that 


(OeffU) = 

where 


4( g 1 u ,.) 2 ( g gx; tt 'l; + a xSxS m k + 2 a xtxf w k + 2 a xixt w xi m x} + 2 a x?xt w XS w xt^ 

(29) 


v- 1 1 

Wi = o w x°! + n w x ° 2 + w xf + L 


(30) 


In the limit of decoupled gauginos, w = w^o ~ w^o fa w ~±, so we obtain [7 


x 2 xf 
2 


(cr e ftv) = 


a HH W 

(^ + l) 2 ’ 


where 


a iiH ~ 


81 g* + 12 g ,2 g 2 + 43// 


„/4 


20487t/i 


w = 


eff 


/^eff \ 

Ms/ 


3/2 


exp 


- x (tiL - 1 
xMs 


(31) 


(32) 


In this case, the relic density can be determined from Eqs. (26) and (31) and the parameter 
space for the DM mass and the mass splitting between Higgsinos and DM is shown in Fig. [3} 
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Below the red region, the relic density is smaller than the lower end of the Planck 3cr values 
so the parameter space in that region is consistent with those obtained for explaining the 
X-ray line in Fig. [2] 

If the Higgsino mass splitting is not ignorable, the contribution of the lighter neu¬ 
tral/charged Higgsino to the effective annihilation cross section gets larger, while the one 
of the heavier Higgsino gets smaller. But, overall, the effective annihilation cross section 
would increase due to the lighter charged Higgsino. Therefore, the difference between the 
averaged Higgsino mass and the singlino can be larger than /x ef j — M$ as shown in Fig. [3] 


5 Collider searches 


When the gauginos are heavy enough, the mass splitting between neutral Higgsinos is of 
about keV scale, being as small as the one between singlinos, and the charged Higgsino 
is almost degenerate in mass with the neutral Higgsinos. On the other hand, the dif¬ 
ference between Higgsino and singlino masses should be less than about 10 GeV for the 
co-annihilation with singlino DM. Thus, neutral Higgsinos and charged Higgsino decay 
dominantly into singlinos with three-body modes such as X 2 Xi u Y Xi 2 X 34 Z* and 
xf as well as the modes containing Higgs fields. Since the singlino Yukawa 

couplings are small for the X-ray line, the charged/neutral Higgsinos decay modes into 
singlino plus off-shcll W/Z leave displaced vertex plus missing energy |7|. When the decay 
length of the charged Higgsino is between about 1 cm and 100 m, there is a bound on the 
mass of the charged Higgsino from the disappearing tracks at the LHC [l6|. In our case, 
when the charged Higgsino decays mostly into singlino plus off-shell W, the decay length 
of the charged Higgsino is about O(mm-rn) as in Ref. |7|. Thus, a certain parameter space 
of a small mass splitting can be constrained. However, the LEP, Tevatron and the LHC 
Run I 17 are not sensitive enough to rule out the neutral Higgsinos. 


When Higgsinos have a sizable mass splitting due to the non-decoupling effect of gaug¬ 
inos, the heavier neutral Higgsino can decay into the lighter neutral or charged Higgsino 
with a sizable branching fraction, and the charged Higgsino can decay into neutral Hig¬ 
gsinos as well. In this case, since Higgsinos have gauge interactions, there is no displaced 
vertex. However, depending on the mass splitting of the Higgsinos, missing energy plus 
collimated leptons at the primary vertices can be a signature. In order for the gauginos 
not to give a large contribution to the singlino mass splitting, their contribution to the 
Higgsino mass splitting is less than 5 GeV. In this case, the situation would be better due 
to larger efficiency of the lepton momentum cuts, as compared to the Higgsinos with keV 
mass splitting. The detailed discussion on the search for almost degenerate Higgsinos is 


outside the scope of this work, so it is left for a future publication 14 . 
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6 Conclusions 


We have considered a Dirac singlet fermion or singlinos with small mass splitting as thermal 
DM in the SUSY Higgs-portal model. In order to explain the X-ray line excess observed 
from the sky, we introduced small singlino Yukawa couplings with Higgses and Higgsinos, 
the SUSY version of Higgs-portal couplings, and showed that the mass splitting of 3.55 keV 
is made and at the same time a tiny magnetic transition dipole moment between the 
Majorana components of the singlino is generated. The singlino mass splitting requires 
gaugino masses to be heavier than about 1 TeV, leading to almost degenerate Higgsinos. 
The thermal production of the singlino DM restricts the Higgsino masses to be not greater 
than about 10 GeV as compared to the singlino masses. New search strategies for almost 
degenerate Higgsinos at the LHC Run II and future colliders are needed to probe the SUSY 
Higgs-portal models. 


Acknowledgments 

The work of HML is supported in part by Basic Science Research Program through the Na¬ 
tional Research Foundation of Korea (NRF) funded by the Ministry of Education, Science 
and Technology (2013R1A1A2007919). MP was supported by World Premier International 
Research Center Initiative (WPI Initiative), MEXT, Japan and acknowledges support from 
the Korea Ministry of Science, ICT and Future Planning, Gyeongsangbuk-Do and Pohang 
City for Independent Junior Research Groups at the Asia Pacific Center for Theoretical 
Physics. 


References 

[1] E. Bulbul, M. Markevitch, A. Foster, R. K. Smith, M. Loewenstein and S. W. Ran¬ 
dall, Astrophys. J. 789, 13 (2014) |arXiv:1402.2301 [astro-ph.CO]]; 

[2] A. Boyarsky, O. Ruchayskiy, D. Iakubovskyi and J. Franse, Phys. Rev. Lett. 113, 
no. 25, 251301 (2014) |arXiv:1402.4119 [astro-ph.CO]]; A. Boyarsky, J. Franse, 
D. Iakubovskyi and O. Ruchayskiy, arXiv:1408.2503 [astro-ph.CO]. 

[3] T. E. Jeltema and S. Profumo, arXiv: 1408.1699 [astro-ph.HE]; A. Boyarsky, 
J. Franse, D. Iakubovskyi and O. Ruchayskiy, arXiv: 1408.4388 [astro-ph.CO]; E. Bul¬ 
bul, M. Markevitch, A. R. Foster, R. K. Smith, M. Loewenstein and S. W. Randall, 
arXiv: 1409.4143 [astro-ph.HE]; T. Jeltema and S. Profumo, arXiv: 1411.1759 [astro- 
ph.HE]; E. Carlson, T. Jeltema and S. Profumo, arXiv: 1411.1758 [astro-ph.HE], 

[4] H. Ishida, K. S. Jeong and F. Takahashi, Phys. Lett. B 732, 196 (2014) 
|arXiv:1402.5837 [hep-ph]]; D. P. Finkbeiner and N. Weiner, arXiv: 1402.6671 [hep- 


12 




ph]; T. Higaki, K. S. Jeong and F. Takahashi, Phys. Lett. B 733, 25 (2014) 
|arXiv:1402.6965 [hep-ph]]; J. Jaeckcl, J. Redondo and A. Ringwald, Phys. Rev. 
D 89, no. 10, 103511 (2014) |arXiv:1402.7335 [hep-ph]]; H. M. Lee, S. C. Park and 
W. I. Park, Enr. Phys. J. C 74, no. 9, 3062 (2014) |arXiv:1403.0865 [astro-ph.CO]]; 

R. Krall, M. Reece and T. Roxlo, JCAP 1409, 007 (2014) |arXiv:1403.1240 [hep-ph]]; 
J. C. Park, S. C. Park and K. Kong, Phys. Lett. B 733, 217 (2014) |arXiv:1403.1536 
[hep-ph]]; M. T. Frandsen, F. Sannino, I. M. Shoemaker and O. Svendsen, JCAP 
1405, 033 (2014) |arXiv:1403.1570 [hep-ph]]; K. Y. Choi and O. Seto, Phys. Lett. B 
735, 92 (2014) |arXiv:1403.1782 [hep-ph]]; S. Baek and H. Okada, arXiv:1403.1710 
[hep-ph]; M. Cicoli, J. P. Conlon, M. C. D. Marsh and M. Rummel, Phys. Rev. D 
90, no. 2, 023540 (2014) [arXiv: 1403.2370 [hep-ph]]; F. Bezrukov and D. Gorbunov, 
Phys. Lett. B 736, 494 (2014) [arXiv: 1403.4638 [hep-ph]]; C. Kolda and J. Unwin, 
Phys. Rev. D 90, no. 2, 023535 (2014) |arXiv:1403.5580 [hep-ph]]; R. Allahverdi, 
B. Dutta and Y. Gao, Phys. Rev. D 89, no. 12, 127305 (2014) |arXiv:1403.5717 
[hep-ph]]; N.-E. Bomark and L. Roszkowski, Phys. Rev. D 90, no. 1, 011701 (2014) 
|arXiv:1403.6503 [hep-ph]]; S. P. Liew, JCAP 1405, 044 (2014) |arXiv:1403.6621 
[hep-ph]]; Z. Kang, P. Ko, T. Li and Y. Liu, arXiv: 1403.7742 [hep-ph]; S. V. Demi¬ 
dov and D. S. Gorbunov, Phys. Rev. D 90, no. 3, 035014 (2014) [arXiv: 1404.1339 
[hep-ph]]; F. S. Queiroz and K. Sinha, Phys. Lett. B 735, 69 (2014) [arXiv: 1404.1400 
[hep-ph]]; E. Dudas, L. Heurtier and Y. Mambrini, Phys. Rev. D 90, no. 3, 035002 
(2014) |arXiv: 1404.1927 [hep-ph]]; K. S. Babu and R. N. Mohapatra, Phys. Rev. 
D 89, no. 11, 115011 (2014) |arXiv:1404.2220 [hep-ph]]; J. M. Cline, Y. Farzan, 
Z. Liu, G. D. Moore and W. Xue, Phys. Rev. D 89, no. 12, 121302 (2014) 
|arXiv:1404.3729 [hep-ph]]; S. Chakraborty, D. K. Ghosh and S. Roy, JHEP 1410 
(2014) 146 |arXiv:1405.6967 [hep-ph]]; K. Cheung, W. C. Huang and Y. L. S. Tsai, 
arXiv: 1411.2619 [hep-ph], 

[5] H. M. Lee, Phys. Lett. B 738 (2014) 118 j arXiv: 1404.5446 [hep-ph]]. 

[6] K. P. Modak, arXiv: 1404.3676 [hep-ph]; C. Q. Geng, D. Huang and L. H. Tsai, 
JHEP 1408 (2014) 086 |arXiv:1406.6481 [hep-ph]]; C. W. Cliiang and T. Yamada, 
JHEP 1409 (2014) 006 j arXiv: 1407.0460 [hep-ph]]; J. M. Cline and A. R. Frey, JCAP 
1410 (2014) 10, 013 |arXiv:1408.0233 [hep-ph]]; S. Baek, arXiv: 1410.1992 [hep-ph]; 

S. Patra, N. Sahoo and N. Sahu, arXiv:1412.4253 [hep-ph]; G. Arcadi, L. Covi 
and F. Dradi; arXiv: 1412.6351 [hep-ph], A. Berlin, A. DiFranzo and D. Hooper, 
arXiv: 1501.03496 [hep-ph], 

[7] A. Falkowski, Y. Hochberg and J. T. Ruderman, JHEP 1411, 140 (2014) 
|arXiv:1409.2872 [hep-ph]]. 

[8] X. Lu, H. Murayama, J. T. Ruderman and K. Tobioka, Phys. Rev. Lett. 112 (2014) 
191803 [arXiv: 1308.0792 [hep-ph]]. 

[9] A. Kaminska, G. G. Ross, K. Schmidt-Hoberg and F. Staub, JHEP 1406 (2014) 153 
|arXiv:1401.1816 [hep-ph]]. 


13 




[10] D. Schmidt, T. Schwetz and T. Toma, Phys. Rev. D 85 (2012) 073009 
|arXiv:1201.0906 [hep-ph]]. 

[11] W. S. Cho, J. H. Huh, I. W. Kim, J. E. Kim and B. Kyae, Phys. Lett. B 687 (2010) 
6 [Erratum-ibid. B 694 (2011) 496] |arXiv:1001.0579 [hep-ph]]. 

[12] H. E. Haber and D. Wyler, Nucl. Phys. B 323 (1989) 267. 

[13] K. Griest and D. Seckel, Phys. Rev. D 43 (1991) 3191. 

[14] H. M. Lee, C. B. Park and M. Park, To appear. 

[15] P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys. 571 (2014) A16 
|arXiv:1303.5076 [astro-ph.CO]]. 

[16] [CMS Collaboration], JHEP 01 (2015) 096 |arXiv: 1411.6006 [hep-ex]]. 

[17] G. Aad et al. [ATLAS Collaboration], Phys. Rev. Lett. 108 (2012) 251801 
|arXiv:1203.1303 [hep-ex]]; G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 721 
(2013) 32 |arXiv:1210.0435 [hep-ex]]; ATLAS collaboration, ATLAS-CONF-2013- 
092, ATLAS-COM-CONF-2013-108; G. Aad et al. [ATLAS Collaboration], JHEP 
1411 (2014) 088 [arXiv: 1409.0746 [hep-ex]]; CMS Collaboration, CMS-PAS-EXO- 
12-037; CMS Collaboration, CMS-PAS-EXO-12-038; CMS Collaboration, CMS-PAS- 
B2G-12-024. 


14 


